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Abstract. We consider the nodal domains of Gaussian random waves in two 
dimensions. We present a method to calculate the distribution of the number of 
nodal domains and the average connectivity with the help of auxiliary Potts-spins. 
An analytical approach could be helpful to decide whether the pattern of nodal 
domains belongs to the universality class of short-ranged percolation. This is not 
completely evident since we find (weak) long-ranged correlations between distant 
avoided nodal intersections. 
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1. Introduction 

In the classical world one can easily distinguish between integrable and chaotic systems 
just by looking at the system's trajectories. Upon increasing h and entering the 
quantum world, the notion of trajectory loses its meaning and more subtle measures 
to characterize quantum chaos are needed. The most common method to reveal 
the quantum chaotic nature of a system is to look at the spectral statistics of its 
Hamiltonian. The spacings of the spectral lines have a Wigner-Dyson distribution 
instead of a Poisson or regular distribution. Recently, Blum et. al. P proposed an 
alternative measure based on the distribution of the number of nodal domains of the 
wave function. To have a specific system in mind we consider a particle in a (two- 
dimensional) quantum billiard. In the absence of a magnetic field one can choose real 
wave functions. The nodal domains are now the regions of constant sign; they are 
bounded by the nodal lines where the wave function vanishes. 

In case of an separable system, one observes a grid of intersecting nodal lines, 
whereas in the chaotic case the nodal domains form a highly irregular pattern and 
nodal intersections are rare [3 01 ■ Similar pattern are also found in Gaussian random 
waves, which are known to be a good model for chaotic wave functions |4j. Here, 
nodal intersections hardly ever occur, i.e. saddle points and nodal lines coincide 
with vanishing probability. Instead, depending on the sign of the wave function right 
at a saddle point, junctions between positive (negative) regions are formed [3 El El 
(avoided intersections), as illustrated in figure ^ Bogomolny and Schmit [3 mimic 
the distribution of nodal domains by the following procedure. They start from a 
checkerboard-like lattice of positive and negative regions, where the (linear) size of 
the squares corresponds to the wavelength of the random waves. At each intersection 
they place with equal probability and independently a junction (see figure ^ which 
either connects the positive or negative regions. Since their model belongs to the 
universality class of critical, short-ranged percolation, they are able to compute the 
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moments of the numbers of nodal domains, the area distribution and corresponding 
(critical) exponents. 

In this article we argue, that the applicability of short-ranged percolation to the 
nodal domains of Gaussian random waves is not entirely obvious. We will demonstrate, 
that distant avoided intersections are in fact correlated. Since the correlations are 
rather weak, they will if at all modify the critical properties of very large systems, 
which cannot be accessed numerically, calling for an analytical scheme to compute the 
number of domains and their (mean) connectivity. We put forward a possible method 
in the second part of the article. We hope that it will ultimately lead to an analytical 
approach to the problem of counting nodal domains. 

2. Long-ranged correlations of avoided intersections in random waves 

Gaussian random waves might be seen as a superposition 4> — J2k '^k4>k of partial 
waves 4>k{'^) with fixed energy k'^ and amplitudes ak with a Gaussian distribution. 
The random waves obey the Helmholtz equation — V^^ = k'^(j) and are characterized 
(in two dimensions) by the correlation function ^ 0] U] 

(0(r)0(O)) = Jo(r) - cos(27rr)/yf for r ^ c» (1) 

where we have set the energy to £^ = fc^ = 1. The Gaussian average is denoted by 
(...). The nodal domains are the (connected) regions of constant sign of the random 
field a{r) — sign(0(r)), the nodal lines are at (f>= 0. 

In the pattern of nodal domains is modelled by a staggered (checkerboard- like) 
lattice of positive and negative regions. With equal probability (and at each node 
independently) either two positive or two negative areas are connected (see figure QJ. 
The equivalent quantity for the random waves is simply the sign of the random field (j) 
right at the saddle points — if the sign of is positive at the saddle point, the positive 
regions are connected and vice versa. 

Therefore, we have to calculate correlations of the density 

X = -sign(0)0(- det(a,9,0)) detid,d,(f>)S\V4>). (2) 




Figure 1. Both types of avoided intersections in the model |7| 
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The (5-function picks up the extrema of cf) (where V0 ~ 0), the determinant of didj4> 
yields the corresponding normahzation, and the theta-function singles out the saddle 
points, which have a negative determinant, see e.g. jH]. 

The calculation of (x(''')x(0)) is done in an expansion in terms of Jo{r). One might 
think of X as a series in <f> and its derivatives (to obtain the series expansion, one has 
to replace at first the 5, 9, and sign- function by their respective smooth counterparts). 
Now we perform the Gaussian average with the help of the Wick theorem. We have 
internal contractions within x(r') and within x(0), and contractions connecting xi'"') 
and x(0)- We sort the expansion by the number of the latter connections. The first 
term of the expansion is obviously 

«.W0,, . /dV /.V ) W-O) (^)+^ . , (3, 

It contains all terms, which decay as Joif) ^ r~^/^, whereas the next term in the 
expansion decays not slower than (Jo(r'))'^ ^ ^.-3/2 ^gjj^Qg jg antisymmetric in the 
term c>c (Jo(r))^ vanishes). 

To obtain the functional derivative of x> we add to an infinitesimal (fixed) 
variation ip ^ (p + rj and calculate (x([</' + ^]) ■'")) to linear order in rj. Since dcf) is 
independent from </> and from dd4> as well, we have 

{X{[<l> + V],r)) - (sign((/. + 77)F(-det(a,9,(/. + 9,9,77))) (<52(V0 + V77)) , (4) 

where F{x) = x9{x). Due to its antisymmetry with respect to (/>—!■ —0 the left 
expectation value vanishes for 77 = 0. Therefore, using the derivatives sign'(a;) = 25{x) 
and F'{x) = e{x) 

Sr^ixir)) = (^2(V</))) 277 ((5(0)F(-detaa0)) 

- (<52(V0)) (sign(0)0(- detddcl))d,dk(l)) t^tkid.dni (5) 

where the calculation of equation (|SJ| is done in [Appendix A| Utilizing once more 
the wave equation Jo = — Jq we obtain finally the large distance behaviour of the 
correlation function 

(x(r)x(O)) = ^Mr) + 0(J^(r)) ^ cos(27rr)/V^. (7) 

The correlations of x are in fact long ranged. The prefactor, however, is extremely 
small, if compared with a naive guess for the x correlation function, namely 

2 2 4 

(sign((/)(r))sign((/)(0))) x (density of saddle points) ~ —Jo{r) x - + . . . (8) 

TT 3 

The prefactor of the naive approximation is 1927r^ « 1895 times larger than the 
prefactor of the exact asymptotic result. The numerical smallness of the x correlations 
is the key to an understanding of the apparent agreement of the results of [7] with the 
real random waves One of the two following scenarios might be realized: Either 
the long-ranged correlations are irrelevant due to their oscillatory character and the 
system belongs to the short-range percolation universality class. Or, the behaviour of 
the system eventually crosses over to a behaviour, which is different from short-range 
percolation. Due to the smallness of iQ, the crossover will take place at very large 
systems (deep in the semiclassical regime) — for small systems, one will see percolative 
behaviour in agreement with [7j. 
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3. Counting nodal domains with the help of Potts-spins 

In this section we study the percolative structure of the nodal domains of random 
waves. To be more specific, we attempt to calculate the moments of the number of 
nodal domains. As a byproduct we will obtain the probability, that two points at 
a distance r belong to the same nodal domain. To achieve this we make use of a 
Potts-model "on top" of the nodal structure. For regularization purposes, the nodal 
domains are placed on a quadratic lattice with a lattice parameter a <C 27r/fc. To each 
node Ti (where i enumerates the nodes) of the lattice we assign beside the sign of the 
random wave ai = sign(0(ri)) a Potts spin Si which takes values 1 . . .q, where q is a 
given natural number. We denote the pairs of adjacent nodes by < i, j >, i.e. X)<i j> 
is the sum over all bonds of the square lattice. We impose now a crucial restriction on 
the Potts-spins: All spins Si, which belong to the same nodal domain, must have the 
same value: Si = Sj for all i, j of the same nodal domain. Spins from different domains 
might have different values. The constraint is implemented in a local manner: The 
product over all lattice bonds 

n C-P" + ^s.J = n (i - - ^-^.)) (9) 

<i,j> ^ ^ <ij> ^ 

is one, if and only if all factors are equal to one. A certain factor is one, if either 
Ui Gj belong to different domains) or if i,j are in the same domain ai = <Tj and 
the corresponding Potts spins are equal 5s.,sj — 1- Let us now sum over all possible 
spin configurations {si}. For each domain, the spins acquire the q different values 
1 ... 5, i.e. each nodal domain contributes a factor of q to the sum 

^[w]=E n (i-^^(i-'5...,))=9^ (10) 

where C is the number of nodal domains. 

Now we perform the average {■ ■ ■) ^ over the Gaussian random field with 
correlation function and obtain the partition function 

^-E(-pf-/3E ^^(i-^-.)l) -(O, (11) 

where we have in addition exponentiated the statistical weig ht ®. Equation (UTI 
might be seen as a kind of (annealed) diluted Potts model similar to the one presented 
in |51 with broken bonds along the nodal lines. Now we take the somewhat formal 
limit |1(J) q — > 1 and find for the mean number of nodal domains (the higher moments 
can be obtained analogously) 



for /3 ^ c». (12) 



The connectivity of the domains can also be obtained with the help of the above 
partition function. Spins on different domains are independent, i.e. the correlation of 
the order parameter 

=^^.,1-1/9 (13) 

vanishes: 

E E('^^-i- V9)('5s„i-i/'z) = o. (14) 

Sk = \...q Si 
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If both spins are constrained to be equal, the correlation between ipi^ and ipi is in fact 
nonzero 

l/qf^l-l/q. 



s 

The correlation function reads now 

E^fc^' n (i- 

{si} <ij> ^ 

ii k,l are disconnected 



(15) 



i±^(l-5. 



q'^ X 



1 — 1/g if /c, ^ are connected 



Now we take the limit q ^ 1 



d 



(q^-q^-')^l, 



9=1 



perform the average (. . over the Gaussian random field 
probability pki that the points k and I are connected 



(16) 
(17) 

and obtain finally the 



Pki 



d_ 

dq 



1 + (TiCTj 



(1 



(18) 



(Ev'fcV'; n 1 

9=1 \{si} <L3> ^ 

The model (jlOII is obviously a direct descendant of the Potts model, interpreted as 
the random-cluster model JO]. In the usual Potts model 

n (1+"'^^..^.) (19) 

{s,} <i,j> 

the two processes of cluster generation and counting are interwoven — the expansion of 
the product generates the different clusters and the summation over the Potts-spins 
counts the number of clusters as in equation (|10|l . In our formula, we average over 
the random waves to generate the nodal domains, and average separately over the 
auxiliary Potts-spins to count the domains. 



4. Conclusion 



We could show, that the avoided intersections of Gaussian random waves show long- 
ranged correlations, which are not taken into account in the simple percolation model 
[7]. It is therefore not evident, whether both models belong to the same universality 
class, i.e. share the same critical exponents. 

In order to investigate the percolative structure of the nodal domains of random 
waves, we have proposed a model, which utilizes Potts-spins to count the number of 
domains. One might easily find variants of this model, e.g. by placing the auxiliary 
spins on the dual lattice in order to count the nodal lines or by using a "height-variable" 
as in ^Jj to obtain the area distribution of the nodal domains. All models suffer from 
the fact, that the auxiliary "spins" have infinite couplings (/3 — > oo in equation l|ll|) '). 
Otherwise the interaction range of the auxiliary variables would be finite, making it 
impossible to decide whether distant regions are connected or not. Unfortunately, 
the infinite couplings complicate calculations and approximations considerably. The 
evaluation of the counting formula is therefore left to future investigations. A helpful 
first step would be the formulation of a continuous field theory corresponding to 
equation Hll|l . Then, the relevance (or irrelevance) of the long-ranged correlations 
((TJ would become manifest. 
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Appendix A. Calculation of the averages in equation (jS)) 

At first we calculate 

- (sign(0)0(- dei^^(j))^^^k4>) = - {5^k/2) {sign{cp)V^ 4> B {{d^dy^f - dMdydy<t>)) 

= (,5,fe/2) (101 e - dA<Pdydycf)) (A.i) 

due to the 0(2) symmetry and due to the fact, that (j) obeys the wave equation —\7'^(f) = 
(j). We introduce the Gaussian variables X = y/%dxdy4>, Y = -J2{dxdx — dydy)<l) and 
Z = {dxdx + dydy)(j) = —0, which are mutually independent random variables with 
second moments (X^) = (Y^) = (Z^) = 1. Then -detddc/) = {X^ +Y^ - 2Z^)/8 
and 

- (sign(0)0(- det 

= {S,k/2){\Z\e{X'+Y'-2Z')) 
A r 

' p2dpsin(z9)di?p|cos(i?)|6'(sin2(z9) -2cos2(79))exp(-pV2) 
dz9 sin(t?)| cos(i9)| 6 (l - 3cos2(i9)) 



2V27r 



V2ttJo 

= ^ f\t\t\e{i-3t')^:^, (A.2) 

where we have introduced spherical coordinates X = /9sini?cos(/3, Y = psin'd simp, 
and Z = pcosi?. The other expectation value is 

2 {S{(p)F{- dot ddcl))) = 2 {5{Z)F {{X^ +Y^ - 2Z^)/8)) 

= (5(Z))((x2+y2))/4= 1 . (A.3) 



With {5'^{V4))) = I/tt we obtain finally 

{x{r)) = (27r)-3/2,^ + [in^j^)-^ V^rj. (A.4) 
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